On the role of initial conditions and final state interactions in ultrarelativistic heavy ion 

collisions 
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We investigate the rapidity dependence of the elliptical flow in heavy ion collisions at 200 GeV 
(cms), by employing a three-dimensional hydrodynamic evolution, based on different initial con- 
ditions, and different freeze-out scenarios. It will be shown that the form of pseudo-rapidity (rj) 
dependence of the elliptical flow is almost identical to space-time-rapidity (t] b ) dependence of the 
initial energy distribution, independent of the freeze-out prescriptions. 



A hydrodynamical treatment of ultrarelativistic 
heavy ion collisions requires thermalized matter as 
an "initial condition" at some early time being of the 
order of a fm/c. Such an initial condition is difficult 
to access theoretically. There are also considerable 
uncertainties concerning the evolution of the system, 
concerning transport coefficients and the equation of 
state. Finally towards the end of the evolutions, it 
seems more and more clear that the system will not 
stay in thermal equilibrium, but interact neverthe- 
less via hadronic rescatterings, before freezing out. 

It is therefore desirable to disentangle the different 
phases, try to understand which kind of observables 
are sensitive to what feature of the model descrip- 
tion. In this paper we are going to investigate the 
role of the initial condition, and we will in particular 
focus on the rapidity dependence. 

We will compare two options for initial condi- 
tions: a parameterization, with parameters chosen 
in order to optimize final results referred to as 
"PAR" thoughout this paper, and an initial condi- 
tion obtained from microscopic approach "EPOS", 
based on the hypothesis that thermalization hap- 
pens very quickly and is achieved at some To. For 
both options, we will perform three-dimensional hy- 
drodynamic calculations, using the same equation of 
state, see For either of the two scenarios, we will 
investigate different freeze-out (FO) scenarios, to be 
discussed later. This modular structure allows us to 
separate initial and final state effects. 

The PAR initial condition has been employed in 
many publications, for details see In case of 
EPOS, the initial scatterings lead to the formation 
of strings, which break into segments, which are usu- 
ally identified with hadrons. When it comes to heavy 
ion collisions, the procedure is modified: one consid- 
ers the situation at an early proper time t , long 
before the hadrons are formed: one distinguishes 



between string segments in dense areas (more than 
some critical density po segments per unit volume), 
from those in low density areas. The high density 
areas are referred to as core, the low density areas 
as corona Q. The corona is important for certain 
aspects, not the ones looked at in this paper. So 
here we simply consider the core part. In any case, 
it is important to note that initial conditions from 
EPOS are based on strings, not on partons. Based 
on the four-momenta of the string segments which 
constitute the core, we compute the energy density 
e(tq,x) and the flow velocity v(to,x). 

Having fixed the initial conditions, the system 
evolves according the equations of ideal hydrody- 
namics, for details see [lj. To have a flexible and 
mudular structure, we make first FO tables (storing 
FO surface and flows) based on hydro calculations 
with PAR and EPOS initial conditions (for given 
Tpo)- We then generate particles EbE from the 
core, using FO tables, based on the Cooper-Frye pre- 
scription. 

In the following figures, "EPOS" refers to the hy- 
drodynamic evolution based on EPOS initial condi- 
tions, "PAR" refers the parameterized initial condi- 
tions of Both calculation use to = 0.6fm/c, and 
the same equation of state. In fig. [IJ we show the 
space-time evolutions of the energy density s, for a 
central Au+Au collision (b = 2.2 fm), for PAR (up- 
per plots) and EPOS (lower plots) initial conditions. 
In each plot, the different curves refer to different 
times, from top to bottom: r = tq, t = tq + At, 
t = to + 2At, etc, with At = 3.9 fm (4 fm) in case 
of PAR (EPOS). The left plots show £ as a function 
of the space-time rapidity 77, for a transverse distance 
r = 0, the middle plots show e as a function of the 
transverse distance r, for 77 = 0, the right ones for 
rj = 2. Full (dotted) curves refer to <p = (<£> = 7r/2). 
The most striking difference between the PAR and 
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Figure 1: Space-time evolutions of the energy density e, for a central Au+Au collision (6 = 2.2 fm), for PAR (upper 
plots) and EPOS (lower plots) initial conditions. In each plot, the different curves refer to different times, from top 
to bottom: r = r , t = r + At, t = r + 2Ar, etc, with At = 3.9 fm (4 fm) in case of PAR (EPOS). The left plots 
show e as a function of the space-time rapidity, for a transverse distance r = 0, the middle plots show e as a function 
of the transverse distance r, for rj = 0, the right ones for r] = 2. Full (dotted) curves refer to if = (ip = 7r/2). 



EPOS curves is the space-time rapidity dependence 
of the energy density, at initial time To (upper curves 
on the left plots) : Whereas the PAR curve is flat up 
to r) = 2, the EPOS curve drops considerably. 

In fig. [21 we show the space-time evolutions of 
the energy density e, for a peripheral Au+Au colli- 
sion (6 = 8.2 fm), for PAR (upper plots) and EPOS 
(lower plots) initial conditions. Again, the most 
striking difference is the stronger space-time rapidity 
dependence of the initial distribution. In transverse 
direction the curves are quite similar, however, the 
EPOS results have a bigger eccentricity, as can be 
seen from the bigger difference between the ip = 
and the tp = ir curves. 

In the following, we will discuss results concerning 
the elliptical flow V2 as a function of the pseudora- 
pidity 77, for Au+Au collisions at 200 GeV. For both 



PAR and EPOS, we show always three curves: 

• a dashed one, representing a pure hydrody- 
namic evolution, with freeze out 100 MeV 
(FO_100); 

• a dotted curve, referring to a pure hydrody- 
namical evolutions, with freeze out at 169 MeV 
(FO_169); 

• a full curve, referring to a hydrodynamical 
evolution with freeze out at 169 MeV, and 
subsequent hadronic cascade, using UrQMD 
(FO_169+cascade). 

We should keep in mind that the critical tempera- 
ture is T c = 170 MeV . In figOU we show first the v 2 
results for the PAR initial conditions, for minimum 
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Figure 2: Space-time evolutions of the energy density e, for a peripheral Au+Au collision (b — 8.2 fm), for PAR 
(upper plots) and EPOS (lower plots) initial conditions. In each plot, the different curves refer to different times, 
from top to bottom: r = to, t — r + At, t — r + 2Ar, etc, with At = 3.9 fm (4 fm) in case of PAR (EPOS). 
The left plots show £ as a function of the space-time rapidity, for a transverse distance r = 0, the middle plots show 
e as a function of the transverse distance r, for rj — 0, the right ones for r) = 2. Full (dotted) curves refer to if = 
(¥> = 7r/2). 



bias (MB) events, as well as for different central- 
ity classes: 3-15%, 15-25%, 25-50%, and 0-40% of 
the most central collisions, compared to data [jj. 
The data show for all centralities a more or less 
pronounced triangular shape, with larger values for 
more peripheral collisions. The calculations show 
as well a very similar shape for the different cen- 
trality classes, for each of the three freeze out op- 
tions. The FO_100 option gives the largest v 2 val- 
ues, the shape is quite flat. The FO_169 option 
(early freeze out) gives significantly smaller i>2 val- 
ues, and the distributions get narrower. Consider- 
ing finally early freeze out with subsequent hadronic 
cascade (FO_169+cascade), leads again to some in- 
crease of V2, but it remains considerably lower than 
the pure hydro results FO_100. 



The theoretical curves are identical to those shown 
in [l[, in case of the FO_100 and FO_169 options. 
Concerning the hadronic rescattering, two different 
approaches have been employed: the UrQMD model 
in this work, and JAM in ref. [l|. The results are 
quite close, though not identical. The JAM curves 
show a slight dip at midrapidity, which is absent 
when using UrQMD. It is difficult to pin down the 
origin of these differences, but it is encouraging that 
the differences are relatively small, and they may be 
considered as the "systematic error" of the theoreti- 
cal treatment of this part. 

In fig. IH we show the corresponding t>2 results 
for the EPOS initial conditions. Here, the shape 
is completely different compared to the PAR initial 
conditions, it is more triangular. The different freeze 
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Figure 3: The elliptical flow vi as a function of the pseudorapidity rj, for the PAR initial conditions, for minimum 
bias (MB) events, as well as for different centrality classes: 3-15%, 15-25%, 25-50%, and 0-40% of the most central 
collisions, compared to data The three curves in each plot refere to different freeze out prescriptions: FO_100 
(dashed), FO_169 (dotted), FO_169+cascade (full). 



out options differ in magnitude (in the same fashion 
as for the PAR results), but they all show a similar 
overall shape. 

How can we understand this large difference in the 
shape of the 77 dependence of V2, between PAR and 
EPOS initial conditions, which is even independent 
of the freeze out prescriptions? From figs. 1 1121 we 
know, that the space time evolution of the energy 
density in PAR and EPOS are quite similar, apart 
of the fact that initially the space-time rapidity ri s 
dependence in PAR is flat, wheras in EPOS it drops 
strongly with rj s . Using a linear scale, initial rj s de- 
pendence is already almost triangular. So there is a 
strong correlation between the r? s width of the initial 
energy density and the pseudorapidity distribution 
of V2 , which is compatible with earlier findings 0| • 

The triangular-like shape of the initial r? s distri- 
bution in EPOS is, on the other hand, partly due to 



the fact that we use strings as a basis of the calcu- 
lation of the initial energy density. Strings always 
strech over a certain range in r? s , with fluctuations 
concerning the length of the string, but always cov- 
ering i] s — 0, leading thus to a triangular-like shape. 

We have seen significant differences in the rapidity 
dependence of the elliptical flow, when choosing dif- 
ferent initial conditions. Do similar differences show 
up in simple rapidity spectra as well? In fig. [5j 
we show the rapidity distributions of pions in cen- 
tral Au+Au collisions, for the PAR initial conditions 
(upper panel), and EPOS initial conditions (lower 
panel). Again we consider the three different freeze 
out scenarios, but they provide quite similar results, 
although the late freeze out (FO_100) gives slightly 
more particles than the two other options, the lat- 
ter ones being almost identical. Comparing PAR 
and EPOS initial conditions, the former one gives a 
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Figure 4: The elliptical flow i>2 as a function of the pseudorapidity 77, for the EPOS initial conditions, for minimum 
bias (MB) events, as well as for different centrality classes: 3-15%, 15-25%, 25-50%, and 0-40% of the most central 
collisions, compared to data Q. The three curves in each plot refere to different freeze out prescriptions: FO_100 
(dashed), FO_169 (dotted), FO_169+cascade (full). 



broader distribution, as expected. But, the differ- 
ence is quite small, much smaller than the difference 
concerning the v2 rapidity dependencies. 

To summarize: we compared two options for ini- 
tial conditions for 3D hydrodynamical calculations 
of AuAu collisions: a parameterization, with param- 
eters chosen in order to optimize final results (PAR) , 
and an initial condition obtained from a microscopic 
approach (EPOS). For both options, we performed 
the hydrodynamic calculations using the same equa- 
tion of state. For either of the two scenarios, we in- 
vestigated three different freeze-out options: a pure 
hydrodynamic evolution, with freeze out 100 MeV, 
a pure hydrodynamical evolution, with freeze out 
at 169 MeV, and a hydrodynamical evolution with 
freeze out at 169 MeV, with subsequent hadronic cas- 
cade. We found a fundamentally different shape of 
the pseudorapidity dependence of W2for the two dif- 



ferent initial conditions PAR and EPOS, indepen- 
dent of the freeze out options. The characteristic 
triangular shape for the V2 results for EPOS initial 
conditions is due to a triangular space-time rapidity 
dependence of the initial energy density. The latter 
fact is mainly due to the fact that strings are used 
to compute initial conditions, not partons. It is in- 
teresting to note that rapidity spectra are much less 
affected by the the choice of initial conditions than 
the V2 rapidity dependence. 
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Figure 5: Rapidity distributions of pions in central Au+Au collisions, for the PAR initial conditions (upper panel), 
and EPOS initial conditions (lower panel). The three curves in each plot refere to different freeze out prescriptions: 
FO_100 (dashed), FO_169 (dotted), FO_169+cascade (full). 
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